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PolyChain 


A  Computer  Program  for  Reliability  Evaluation  of 
Large-Scale  Undirected  Networks  via  Polygon-to- 

Chain  Reductions 


by 

Nauricio  G.  C.  Resende 

Operations  Research  Center 
University  of  California,  Berkeley 


1 .  Introduction 


The  evaluation  of  network  reliability  [A83],  is  of  great 
importance  in  engineering.  It  is  both  an  important  design 
and  operations  parameter  in  such  systems  as  communications, 
power,  and  computer,  to  name  a  few. 

Satyanar ayana  and  Wood  [S82],  [W82]  introduced  a  linear  time 
algorithm  for  computing  the  reliability  of  a  network  with  an 
underlying  series-parallel  structure.  This  algorithm 
employs  the  new  polygon-to-chaln  reductions .  They  also 
present  an  extension  to  the  algorithm  which  enables  the  gen¬ 
eration  of  a  reduced  network  when  the  input  network  is  not 
totally  reducible. 

This  report  discusses  the  design  and  implementation  of 
PolyChain  (Polygon-to-Chain ) ,  a  portable  FORTRAN  program  for 
K-terminal  network  reliability  evaluation  via  polygon-to- 
chain  reductions.  Doubly  linked  multilists  and  linear 
stacks  are  used  in  this  implementation  of  the  algorithm  in 
such  a  way  that  large-scale  undirected  networks  are  treated 
efficiently.  PolyChain  is  modular,  thoroughly  commented, 
has  input  data  consistency  tests,  and  is  structured  whenever 
possible,  to  facilitate  maintenance.  There  are  four  output 
options,  which  produce  flexible,  informative  reports. 

Section  2  briefly  discusses  some  of  the  theoretical  results 
of  polygon-to-chain  reductions.  In  section  3*  a  system's 
manual  describes  the  algorithm's  implementation  in  FORTRAN. 
Section  4  consists  of  a  user's  guide  to  PolyChain.  In  sec¬ 
tion  5,  a  number  of  large-scale  problems  problems  are  tested 
on  PolyChain,  illustrating  the  code's  performance  capabili¬ 
ties.  Conclusions  and  recommendations  are  made  in  section 


s  s 


3 


2.  Polygon-to-Chaln  Reductions  and  Series-Parallel  Graphs 


For  a  complete  discussion  of  polygon-to-chain  reductions  and 
series-parallel  graphs,  see  [S82],  C  W  8  2  3  -  In  this  section, 
the  K-teroinal  network  reliability  problem  will  be  defined, 
dasic  definitions  will  be  made  so  that  the  algorithm  of 
Satyanarayana  and  Wood  can  be  presented. 

2.K  K-Terminal  Reliability 

Let  Gs(V,E)  be  a  graph,  where  V  and  E  are  respectively  G's 
vertex  set  and  edge  set.  We  assume  that  vertices  function 
perfectly  and  that  edges  have  a  probability  of  functioning, 
which  may  be  less  than  1.  Edge  e(i>  has  probability  p(i)  of 
functioning  and  q(i)  =  1  -  p(i)  of  not  functioning.  Associ¬ 
ated  with  G  is  a  subset  of  V  denoted  by  K.  These  dis¬ 
tinguished  vertices  are  called  the  K-vertices  of  G.  G(K)  is 
graph  G  with  K  specified.  The  K-terminal  reliability  of 
G(K),  R[G(K)3,  is  the  probability  that,  at  a  given  time  t, 
all  K-vertices  are  connected  by  working  edges. 

2.2.  Simple  Reductions 

To  compute  the  K-terminal  reliability  of  G,  it  is  desirable 
to  reduce  the  size  of  G.  There  are  three  well  known  simple 
reductions:  parallel  reduction,  series  reduction  and 
degree-2  reduction.  In  parallel  reduction ,  parallel  edges 
e(a)s(u,v)  and  e(b)=(u,v)  are  replaced  by  edge  e(c)=(u,v) 
with  edge  probability  p(c)  =  1  -  q(a)qCb).  In  series  reduc¬ 
tion  ,  edges  e(a)=(u,v)  and  e(b)s(v,w)  are  replaced  by  edge 
e(c)=(u,w)  with  probability  p(c)  =  pCa)p(b).  Let  u,v,  and  w 
be  edges  in  set  K.  Let  deg(v)  =  2  and  consider  two  edges 
e(a)=(u,v)  and  e(b)=(v,w),  such  that  u  is  not  equal  to  w.  A 
degree-2  reduction  substitutes  edges  e(a)  and  e(b)  by  edge 
e(c)s(u,w)  with  p( c) =p( a) p ( b ) /[ 1 -q ( a)q( b) ]  and  R[G(K)]=[1- 
q( a)q(b) ]R [G (K-v) ] ,  where  G(K-v)  is  the  new  reduced  graph. 
These  simple  reductions  are  examples  of  reliability¬ 
preserving  reductions  ,  where  a  graph  G  is  replaced  by  a 
reduced  graph  G',  and  R (G )  =  OMEGA*R (G  ’ ) ,  where  OMEGA  depends 
on  the  reduction. 

£. 3 .  Series-Parallel  Reducible  and  Complex  Graphs 

A  series-parallel  graph  is  a  graph  that  can  be  reduced  to  a 
tree  after  successive  series  and  parallel  reductions.  For  a 
given  series-parallel  graph  G,  G(K)  may  or  may  not  be 
reduced  to  a  single  edge  by  successive  simple  reductions. 
This  will  depend  on  the  elements  of  set  K.  G(K)  is  series- 
parallel  reducible  if  it  can  be  reduced  to  a  single  edge  by 
successive  simple  reductions.  It  is  series-parallel  complex 
if  it  is  not  series-parallel  reducible. 


2 . 4_.  Chains  and  Polygons 

A  chain  X  is  an  alternating  sequence  of  distinct  vertices 
and  edges,  such  that  all  vertices,  except  for  the  two  end¬ 
points,  have  degree  2.  Let  XI  and  X2  be  two  distinct  chains 
with  common  endpoints.  The  union  of  XI  and  X2  constitutes  a 
polygon . 


2.5. 


ton-to-Chain  Reductions 


The  main  result  in  [S82],  [W82]  follows.  Let  G(K)  be  a 
graph  that  admits  no  simple  reductions.  If  G(K)  contains  a 
polygon,  it  is  one  of  the  eight  types.  A  reliability 
preserving  reduction  permits  the  replacement  of  the  polygons 
by  chains  with  the  new  edge  reliability  computed  according 
to  given  formulas. 

2.6.  An  0 ( | E | )  Algorithm  for  Series-Parallel  Complex  Graphs 


An  0(|Ei)  algorithm  for  a  nonseparable  series-parallel  com¬ 
plex  graph  is  presented  in  [S82],  [W82].  PolyChain  is  a 
direct  implementation  of  that  algorithm  utilizing  extension 
4'  in  place  of  step  4. 


2*  System’s  Manual 


In  this  section,  the  code  is  briefly  described.  First,  pro¬ 
gramming  is  commented  and  data  structures  discussed.  Next, 
data  consistency  tests  are  considered.  A  data  dictionary  is 
presented.  The  program's  COMMON  blocks  and  all  procedures 
are  described,  along  with  their  corresponding  inputs  and 
outputs.  Finally,  we  comment  on  the  code's  debugging 
feature . 

3..JL*  Program  Design  and  Implementation 


Programming 

PolyChain  is  designed  to  be  portable.  All  algorithmic  rou¬ 
tines  are  written  in  standard  FORTRAN  IV.  Only  output 
related  code  is  system  dependent.  The  code  is  intended  to 
be  friendly  to  both  the  person  using  it  in  an  application 
and  the  one  extending  or  maintaining  it.  The  program  is 
modular,  has  format  free  input  and  informative  outputs. 

3..1..2.  Data  Structures 

The  algorithm  of  Satyanarayana  and  Wood  manipulates 
undirected  networks.  These  networks  can  be  represented  in 
many  ways  in  a  computer  code.  Matrix  representation  has 
many  drawbacks.  These  matrices  are,  in  practice,  extremely 
sparse.  Densities  of  less  than  one  percent  are  not  uncom¬ 
mon.  Matrix  representation  is  inefficient  both  with  respect 
to  core  usage  and  execution  times.  For  a  static  network,  an 
efficient  representation  is  a  packed  matrix.  In  the  algo¬ 
rithm  considered  here,  the  network  is  dynamic.  Helgason  and 
Kennington  [H80],  Thesen  [T78],  Berztiss  [B75],  among  many, 
discuss  efficient  network  representations  using  linked  list 
data  structures.  Thesen  [T78]  and  Berztiss  [B751  discuss 
the  implementation  of  linked  lists  in  FORTRAN. 

During  the  algorithm's  reduction  process,  one  or  more  edges 
or  vertices  are  removed  from  the  network.  In  the  data 
structure,  this  corresponds  to  removing  elements  from  the 
lists.  This  process  is  repeated  frequently  in  the  algo¬ 
rithm.  Doubly  linked  lists  [K733,tB753,  require  more  core, 
but  are  more  efficient  when  many  element  deletions  are 
required.  They  are  implemented  in  PolyChain.  Each  vertex 
has  a  list  of  adjacent  vertices,  which  besides  indicating 
whio  '  vertices  are  adjacent  to  it,  also  provides  information 
det-  mining  whether  the  vertex  and  its  adjacent  vertices 
be  jng  to  set  K.  For  every  element  of  the  list,  there  is  a 
pointer  indicating  the  address  where  information  about  the 
edge,  corresponding  to  these  two  vertices,  is  kept.  Figure 
I  illustrates  this  multilist  data  structure  for  a  small 


network 


3,.1_.2.1_.  Data  Structure  Implementation 

Next,  we  describe  the  FORTRAN  arrays  used  to  implement  the 
code's  data  structures. 

PTRADJ(i)  points  to  the  beginning  of  the  list  of  vertices 
adjacent  to  vertex  i.  If  PTRADJ(i)  is  positive,  vertex  i 
does  not  belong  in  set  K.  If  PTRADJ(i)  is  negative,  then 
vertex  i  is  in  set  K.  PTRADJ(i)  may  also  be  null.  This 
indicates  that  vertex  i  has  been  removed  from  the  network. 
Array  ADJVRT(*)  is  the  principal  information  element  in  the 
list  structure.  It  contains  an  adjacent  vertex  to  the  ver¬ 
tex  represented  by  the  list  in  which  it  resides.  If 
ADJVRT(*)  is  positive,  that  vertex  is  not  in  K,  and  if  it  is 
negative,  the  vertex  is  in  set  K.  LNKDWN(i)  points  to  the 
next  element  (downwards)  in  the  list  of  adjacent  vertices. 
If  LNKDWN(*)  =  0,  this  element  is  the  last  element  in  the 

list.  LNKUP(*)  points  to  the  element  just  above  it  in  the 
adjacent  vertices  list.  If  LNKUP(*)  =  0,  then  this  element 
is  the  first  in  the  list.  EDGPRBC*)  contains  the  edge  reli¬ 
ability.  LNKEDG(*)  points  to  the  position  in  array 


EDGPRB(*)  corresponding  to  edge  [ PTRAD J ( i ), AD JVRT (*)] . 
AVSADJ  is  a  pointer  to  the  beginning  of  the  list  of  avail¬ 
able  space  [B75]. 

The  algorithm  also  requires  data  structures  for  both  the  T 
list  and  the  chain.  Since  the  order  of  the  vertices  in  both 
the  T  list  and  the  chain  is  irrelevant  to  the  code,  linear 
stacks  are  used  to  represent  both  structures.  TLIST(*)  con¬ 
tains  the  stack  of  vertices,  which  are  elements  of  the  T 
list.  If  TLIST(*)  is  positive,  then  that  vertex  is  not  in 
set  K.  If  it  is  negative,  the  vertex  is  in  set  K.  HOP 
points  to  the  top  of  the  TLIST(*)  stack.  CHAINC*)  is  a 
stack  of  vertices  belonging  to  a  given  chain.  TOP  points  to 
its  top.  As  with  TLIS T ( * )  ,  if  CHAIN(»)  is  positive,  then 
that  vertex  is  not  in  set  K.  If  it  is  negative,  the  vertex 
is  in  set  K.  ONLIST(i)  =  .TRUE,  implies  that  vertex  i  is  on 
the  T  list,  and  ONLIST(i)  =  .FALSE.  implies  the  contrary. 
DEG(i)  indicates  the  degree  of  vertex  i.  Figure  II  illus¬ 
trates  the  use  of  some  of  these  arrays  in  network  represen¬ 
tation  . 
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Figure  II  -  FORTRAN  Implementation  of  List  Structure 


3.. 2.  Input  Data  Consistency  Tests 

As  described  later,  the  input  of  a  vertex  in  set  K  will 
require  that  its  vertex  number  come  preceded  by  a  negative 
sign.  PolyChain  tests  consistency  of  the  signs  of  the  ver¬ 
tices,  and  in  this  manner  will  cease  execution  if  there  is 
any  inconsistency  in  the  sign  of  a  vertex  number.  Edge 
reliability  is  also  tested  for  values  less  than  zero  or 
greater  than  one. 


^.3,.  Data  Dictionary 

Next,  the  code's  variables  are  listed,  with  a  brief  descrip¬ 
tion  of  each  one. 


AD JVRT ( * ) 

ALFA 

AUX 

AUXPTR 

AUXTOP 

AVSADJ 

BETA 

BOTTOM 

CARDE 

CAR  DK 

CAR  DL 

CARDV 

CCHAIN 

CHAIN  ( * ) 

CTEST 

D 1 

D2 

DATE 

DEG (* ) 

DELTA 

DENDEN 

DENSTY 

ECHOIN 

ECHOUT 

EDGOUT 

EDGPRB ( • ) 

EGOUTB 

EGOUTC 

ERRO 

FIRST 

FOUND 

HOUR 

GAMMA 

IN 

INON 
IOUT 
IPTR 
IPTR  1 


vertex  adjacent  to  vertex  whose  list  it  is  on 
reduction  parameter 
auxiliary  variable 
auxiliary  pointer 

pointer  to  next  to  the  la3t  element  of  a  chain 

pointer  to  beginning  of  list  of  available  space 

reduction  parameter 

pointer  to  bottom  of  list 

cardinality  of  set  E 

cardinality  of  set  K 

size  of  list  after  input 

cardinality  of  set  V 

absolute  value  of  CHAIN 

chain  stack 

consistency  test  variable 

degree  of  vertex  VI 

degree  of  vertex  V2 

character  string  contain!  \  g  date 

degree  of  vertex 

reduction  parameter 

denominator  in  network  density  computation 

input  network  density 

indicator  of  input  echo 

indicator  of  reduced  network  output 

edge  selected  to  leave  network 

edge  reliability 

selected  edge  to  be  deleted  from  network 

selected  edge  to  be  deleted  from  network 

indicator  of  error  in  some  routine 

key  indicating  first  passage  in  output  routine 

indicator  of  whether  or  not  a  polygon  was  found 

character  string  containing  hour 

reduction  parameter 

value  of  FORTRAN  input  file 

key  indicating  input  echo  option 

value  of  FORTRAN  output  file 

pointer 

pointer 
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IPTR2 

IPTRK 

IPTRNK 

LINECT 

LNKC 

LNKDWN  ( * ) 
LNKEDG ( * ) 
LNKUP ( * ) 

H 

HAXCHN 

NAXEDG 

MAXLST 

MAXVRT 

MCARDE 

HCAROK 

MCARDV 

MXSTKT 

NEWEDG 

NR  1 

NR  2 

NR  3 

NR  4 

NR  5 

NR  6 

NR7 

NR8 

NSR 

ND2R 

NUMDEN 

OMEGA 

ONLIST (*) 

OUTON 

OUTPUT 

PA 

PAGE 

PB 

PC 

PD 

PE 

PF 

POINT 

PR 

PROB 

PS 

PT 

PT  1 3 
PT23 
PTR 

PTRADJ ( • ) 

PTRB 

PTRC 

PTRD2 

PTRDG2 

PTRK 


pointer 

pointer  to  vertex  in  set  K 
pointer  to  vertex  not  in  set  K 
output  line  counter 
pointer 

pointer  to  next  element  on  list 
pointer  to  cor respond ing  edge 
pointer  to  element  above  in  list 
product  of  all  omegas,  see  [W82] 
maximum  number  of  elements  in  chain  stack 
maximum  number  of  edges 

maximum  number  of  elements  in  adjacent  vertices  list 
maximum  number  of  vertices 

cardinality  of  set  E  at  start  of  procedure 

cardinality  of  set  K  at  start  of  procedure 

cardinality  of  set  V  at  start  of  procedure 

maximum  number  of  elements  in  T  list  stack 

pointer  to  new  edge 

counter  of  reductions  of  type  1 

counter  of  reductions  of  type  2 

counter  of  reductions  of  type  3 

counter  of  reductions  of  type  4 

counter  of  reductions  of  type  5 

counter  of  reductions  of  type  6 

counter  of  reductions  of  type  7 

counter  of  reductions  of  type  8 

counter  of  series  reductions 

counter  of  degree  2  reductions 

numerator  in  network  density  computation 

reduction  parameter 

indicator  of  presence  in  T  list 

key  Indicating  reduced  graph  output  on  report 

final  solution 

reliability  of  edge  a 

output  page  counter 

reliability  of  edge  b 

reliability  of  edge  c 

reliability  of  edge  d 

reliability  of  edge  e 

reliability  of  edge  f 

pointer 

updated  edge  reliability 

reliability  of  edge  (VI ,V2) 

updated  edge  reliability 

updated  edge  reliability 

pointer 

pointer 

pointer 

pointer  to  beginning  of  list  of  adjacent  vertices 

pointer 

pointer 

pointer  to  VD2 

pointer  to  VDG2 

pointer  to  vertex  in  set  K 


PTRNK 

PTR  V 

PTRV1 

PTRV2 

PV  1 

PV2 

PTVK 

PTVNK 

QA 

QB 

QC 

QD 

QE 

QF 

TEST 

TIME ( * ) 

TLIS  T ( * ) 

TOP 

TTOP 

TYPE  1 

TYPE2 

USEOGE 

USVRTX 

USLIST 

V 

VI 

V  2 
V3 

V3V1 

V3V2 

V4 

V5V1 

V5V2 

VD2 

VDG  2 

VK 

VNK 

VRTX 

VV 

VV 1 

VV2 

VV  3 

VV3V1 

VVDG2 

VVNK 

W 

X 

XX 

V 

YY 

YEAR 


pointer  to  vertex  not  in  set  K 
pointer 

pointer  to  vertex  adjacent  to  V 

pointer  to  vertex  adjacent  to  V 

pointer 

pointer 

pointer 

pointer 

failure  probability  of  edge  a 

failure  probability  of  edge  b 

failure  probability  of  edge  c 

failure  probability  of  edge  d 

failure  probability  of  edge  e 

failure  probability  of  edge  f 

key  for  debugging  feature 

execution  time 

T  list  stack 

top  of  chain  stack 

top  of  TLIST  stack 

vertex  type  indicator  for  output 

vertex  type  indicator  for  output 

core  usage  for  edge 

core  usage  for  vertex 

core  usage  for  lists 

vertex 

vertex 

vertex 

vertex  adjacent  to  V2 

vertex  adjacent  to  VI 

vertex  adjacent  to  V2 

vertex 

vertex  adjacent  to  vertex  adjacent  to  VI 

vertex  adjacent  to  vertex  adjacent  to  V2 

vertex  with  degree  2 

vertex  with  degree  greater  than  2 

vertex  in  set  K 

vertex  not  in  set  K 

vertex 

absolute  value  of  V 
absolute  value  of  VI 
absolute  value  of  V2 
absolute  value  of  V3 
absolute  value  of  V3V1 
absolute  value  of  VDG2 
absolute  value  of  VNK 
vertex 
vertex 

absolute  value  of  X 
verte  x 

absolute  value  of  Y 

character  string  containing  year 


COMMON  Blocks 
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Most  of  the  variables,  data  structures,  and  parameters  in 
PolyChain  are  shared  by  the  subroutines  through  COMMON 
blocks.  All  of  the  dimensions  of  these  COMMONS  are  deter¬ 
mined  in  SUBROUTINE  1NILST,  through  the  initialization  of 
variables  MAXEDG ,  MAXVRT ,  and  MAXCHN.  We  will  use  these 
values  when  describing  the  blocks.  Below  is  a  list  of  all 
COMMON  blocks  in  the  code. 


COMMON/BLKOI / 
COMMON /BLK02 / 
C0MM0N/BLK21  / 
COMMON /BLK03 / 
COMMON /BLK 04 / 

COMMON /BLK05 / 
COMMON /BLK06/ 
COMMON/BLK07 / 
COMMON /BLK 08/ 
COMMON /BLK09/ 
COMMON /BLK1 0/ 
COMMON /BLK 11/ 
C0MM0N/BLK12/ 
COMMON /BLK 14/ 
COMMON/BLK 15/ 
COMM0N/BLK30/ 
C0MM0N/BLK31 / 
COMMON/BLK32/ 
COMM ON /BLK 33/ 
C0MM0N/BLK31*/ 
COMMON/BLK35 / 


DEG (MAXVRT  ) 

PTRADJ (MAXVRT) .ADJVRT (2»MAXEDG  ) ,AVSADJ 

LNKDWN (2*MAXEDG ) , LNKUP (2*MAXEDG ), LNKEDG ( 2  * MAXEDG ) 

EDGPRB ( MAXEDG ) ,EDGNUM( MAXEDG) 

TLIST (MAXVRT ), CHAIN (MAXCHN  )  ,TTOP,TOP, 

ONLIST (MAXV  RT ) 

MAXEDG, MAXVRT, MAXLST,MXSTKT, MAXCHN 
CARDE , CARDV , CARDK 
M 

IN.IOUT 

FOUND 

IPTR1 ,IPTR2,IPTR,V1 ,V2,D1 ,D2 

PTRV 1 , PTRV2 , PTR V 

V3V1  ,V3V2,  V3,Vi|,V5V1  ,V5V2 

PT 1 3 , PT23 

P315 

NR1 , NR  2, NR  3, NR  4, NR  5, NR  6, NR  7, NR  8, NS R,ND2R 
MCARDE , MCA RDV,M CARDK 
DATE, YEAR, HOUR 
FIRST 

TIME (2 ) .PAGE 
INON , OUTON 


3,.5_.  Description  of  Subroutines 

Next,  we  present  all  SUBROUTINES  in  the  code.  Each  pro¬ 
cedure  will  be  briefly  described  together  with  its  input  and 
output . 


3..5.J..  MAIN  Program 


Description  The  MAIN  program  controls  the  basic  steps  in 
the  algorithm.  From  it,  subroutines, 
representing  the  various  steps  of  the  algo¬ 
rithm,  are  called.  The  algorithm  initiates 
and  terminates  execution  in  the  MAIN  program. 


Input 


An  undirected  network  and  edge  reliabilities. 


Output  Network  reliability,  when  network  is  com¬ 

pletely  reduced.  Reduced  network  and  M,  when 
total  reduction  is  not  possible. 


3.5.2.  SUBROUTINE  INILST 

Description  This  subroutine  initializes  the  code's  most 

important  variables.  It  creates  the  list  of 
available  space  [B75],  which  is  the  first  step 
in  the  initialization  of  the  list  structures. 

Input  Variables  and  arrays  not  yet  initialized. 

Output  Initialized  variables  and  data  structures. 

The  list  of  available  space  is  also  returned, 
ready  to  be  loaded  in  SUBROUTINE  INDATA ( V ) . 

3.5.3.  SUBROUTINE  INDATA(V) 

Description  This  subroutine  reads  the  input  network.  Out¬ 
put  options  are  read  first.  Each  edge  is  read 
on  a  separate  line.  Edges  are  indicated  by 
their  corresponding  vertices  and  reliability. 
A  vertex  number  is  preceded  by  a  minus  sign  if 
the  vertex  is  in  set  K.  Consistency  tests  are 
performed  on  the  input  data,  and  if  any  incon¬ 
sistency  is  reported,  execution  is  terminated. 
Network  multilist  structure  is  loaded.  Cardi¬ 
nalities  and  degrees  are  computed.  Finally, 
the  input  network  is  optionally  printed. 

Input  Free  format  input.  Output  options.  Network 

and  edge  reliabilities.  No  flags  are  needed 
to  indicate  end-of-file. 

Output  The  loaded  multilist  data  structure.  Vertex 

degrees.  The  input  network  is  printed  as  a 
report  if  it  is  so  desired. 

3.5.4.  SUBROUTINE  HEADER 

Description  This  subroutine  prints  the  page  header  of  the 
input  network  output. 

Input  Page  number,  date,  year,  and  hour. 

Output  Page  header  with  incremented  page  number,  and 

date,  year,  and  hour  of  execution. 

3.5.5.  SUBROUTINE  CONLST 

Description  This  subroutine  constructs  the  initial  T  list. 

Input  Degrees  of  all  vertices. 

Output  A  stack,  TLIST(*),  containing  the  vertices  on 

the  T  list.  A  pointer  TTOP  indicating  the  top 
of  stack  TLIST ( * )  .  An  array,  ONLIST(«), 
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indicating  if  a  vertex  is  on  the  T  list. 

3.5.6.  SUBROUTINE  SERIER(V) 

Description  This  subroutine  performs  a  series  reduction  on 
vertex  V  not  in  set  K. 

Input  Vertex  V  and  the  multilist  structure. 

Output  The  updated  multilist  structure,  with  V  and 

both  of  its  edges  deleted,  and  with  a  new  edge 
inserted.  This  new  edge  has  its  reliability 
computed.  New  cardinalities  of  V  and  E. 

3.5.7.  SUBROUTINE  DEG2R ( V  ) 

Description  This  subroutine  performs  a  degree  2  reduction 
on  vertex  V  in  set  K. 

Input  Vertex  V  and  the  multilist  structure. 

Output  The  updated  multilist  structure,  with  V 

deleted,  along  with  both  of  its  edges,  and 
with  a  new  edge  inserted.  This  new  edge  has 
its  reliability  computed.  New  cardinalities 
of  V  and  E.  The  updated  value  of  M. 

3.5.8.  SUBROUTINE  SEARCH(V) 

Description  This  subroutine  controls  the  process  of 

searching  for  a  polygon  emanating  from  vertex 
V,  and,  if  one  is  found,  reducing  it  to  a 
chain . 

Input  Vertex  V  and  the  multilist  structure. 

Output  The  updated  multilist  structure,  with  a 

polygon-to-chain  reduction  performed.  A  chain 
contained  in  stack  CHAIN(*).  Updated  cardi¬ 
nalities  and  degrees.  Updated  N. 

2.5.9.  SUBROUTINE  TYPE8(V) 

Description  This  subroutine  searches  for  a  polygon  of  type 
8  (two  parallel  edges),  emanating  from  vertex 
V.  If  found,  a  reduction  is  performed. 

Vertex  V.  The  multilist  structure. 

The  updated  multilist  structure.  A  chain  con¬ 
sisting  of  two  vertices  and  one  edge.  An 

indicator  (FOUND)  is  activated  if  a  polygon 
was  found. 


Input 

Output 


3.5.10.  SUBROUTINE  DELETE (V  ,  PTR  ) 

Description  This  subroutine  deletes  the  element  pointed  to 
by  PTR  from  vertex  V's  adjacent  vertices  list. 
Three  cases  are  considered.  The  first,  when 
the  element  is  first  in  the  list.  The  second, 
when  it  is  last  in  the  list.  The  last,  when 
the  element  is  in  the  middle  of  the  list.  In 
each  case,  the  element  is  deleted  by  a  dif¬ 
ferent  set  of  of  commands. 

Input  The  multilist  structure.  Vertex  V.  Pointer 

PTR. 

Output  The  updated  multilist  structure  without  the 

specified  element. 

3.5.H.  SUBROUTINE  FINDK  ( V  ) 

Description  This  subroutine  is  called  from  control  routine 
SEARCH(V),  when  vertex  V  is  in  set  K.  It  con¬ 
trols  the  search  for  a  polygon  and,  when  one 
is  found,  its  reduction  to  a  chain.  Only 
polygons  of  types  2  or  3  may  be  found  emanat¬ 
ing  from  a  vertex  V  in  set  K. 

Input  The  multilist  structure.  Vertex  V. 

Output  This  routine  outputs  pairs  of  adjacent  ver¬ 

tices  to  vertex  V,  along  with  their  degrees, 
to  more  specialized  search  routines  (POLY2CV) 
and  P0LY3(V)>.  It  also  outputs  the 
corresponding  pointers  to  these  vertices  in 
vertex  V's  list.  If  a  polygon  is  found,  the 
indicator  FOUND  is  returned  "on"  to  subroutine 
SEARCH  ( V  ) . 

1‘5.12.  SUBROUTINE  P0LY2(V ) 

Description  This  routine  checks  if  vertex  V,  and  the  pair 
of  adjacent  vertices  specified  in  subroutine 
FINDK(V)  form  a  polygon  of  type  2.  If  a 
polygon  of  type  2  is  found,  this  routine  calls 
routine  REDUC2(V),  which  performs  the 
polygon-to-chain  reduction. 

Input  The  multilist  structure.  Vertex  V  and  a  pair 

of  vertices  adjacent  to  V,  along  with  their 
corresponding  degrees.  Pointers  to  the  pair 
of  adjacent  vertices,  in  V's  list. 

Output  Indicator  FOUND  is  activated  if  a  polygon  of 

type  2  is  found. 
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3.5.13.  SUBROUTINE  REDUC2(V ) 

Description  This  routine  is  called  from  subroutine 
P0LY2(V),  when  a  polygon  of  type  2  is  found. 
It  reduces  the  polygon  found  to  a  chain  of 
three  vertices. 

Input  The  multilist  structure.  The  vertices  of  the 

polygon,  and  their  pointers. 

Output  The  updated  multilist  structure.  Updated 

degrees  and  cardinalities.  A  chain  in  stack 
CHAIN ( * ) .  Updated  value  of  M. 

SUBROUTINE  PQLY3(V ) 

Description  This  routine  checks  if  vertex  V,  and  the  pair 

of  adjacent  vertices  specified  in  subroutine 
FINDK(V)  form  a  polygon  of  type  3.  If  a 
polygon  of  type  3  is  found,  this  routine  calls 
routine  REDUC3(V),  which  performs  the 
polygon-to-chain  reduction. 

Input  The  multilist  structure.  Vertex  V  and  a  pair 

of  vertices  adjacent  to  V,  along  with  their 
corresponding  degrees.  Pointers  to  the  pair 
of  adjacent  vertices,  in  V's  list. 

Output  Indicator  FOUND  is  activated  if  a  polygon  of 

type  3  is  found. 

1*5.15.  SUBROUTINE  REDUC3(V ) 

Description  This  routine  is  called  from  subroutine 

P0LY3(V),  when  a  polygon  of  type  3  is  found. 
It  reduces  the  polygon  found  to  a  chain  of 
three  vertices. 

Input  The  multilist  structure.  The  vertices  of  the 

polygon,  and  their  pointers. 

Output  The  updated  multilist  structure.  Updated 

degrees  and  cardinalities.  A  chain  in  stack 
CHAIN(*).  Updated  value  of  M. 

3.5.16.  SUBROUTINE  FINDNK(V) 

Description  This  subroutine  is  called  from  control  routine 
SEARCH(V),  when  vertex  V  is  not  in  set  K.  It 
controls  the  search  for  a  polygon  and  its 
reduction  to  a  chain,  when  a  polygon  is  found. 
If  vertex  V  is  not  in  set  K,  any  type  of 
polygon  can  emanate  from  it. 
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Input  The  multilist  structure.  Vertex  V. 

Output  This  routine  outputs  to  more  specialized 

search  routines  [POLY2P(V),  PY3467(V),  and 
P0LY15(V)]  pairs  of  adjacent  vertices  of  ver¬ 
tex  V,  and  their  degrees.  It  also  outputs  the 
corresponding  pointers  to  these  vertices  in 
vertex  Vs  list.  If  a  polygon  is  found,  the 
indicator  FOUND  is  returned  "on"  to  subroutine 
SEARCH. 

2* 5. 17..  SUBROUTINE  PY3467 (V) 

Description  This  routine  checks  if  vertex  V,  and  the  pair 
of  adjacent  vertices  specified  in  subroutine 
FINDNK(V)  form  a  polygon  of  type  3,4,6,  or  7. 
If  a  polygon  of  type  3,4,6,  or  7  is  found, 
this  routine  calls  routine  REDC3P(V), 
REDUC4 ( V ) ,  REDUC6 ( V  ) ,  or  REDUC7(V),  respec¬ 
tively,  which  performs  the  polygon-to-chain 
reduction . 

Input  The  multilist  structure.  Vertex  V  and  a  pair 

of  vertices  adjacent  to  V.  Pointers  to  the 
pair  of  adjacent  vertices,  in  V's  list. 

Output  Indicator  FOUND  is  activated  if  a  polygon  of 

type  3,4,6,  or  7  is  found.  All  vertices  in 
the  polygon  found,  together  with  pointers 
indicating  their  addresses,  are  passed  to  the 
polygon-to-chain  reduction  routines. 

2.5.18.  SUBROUTINE  REDC3P (V  ) 

Description  This  routine  is  called  from  subroutine 

PY3467(V),  when  a  polygon  of  type  3  is  found. 
It  reduces  the  polygon  found  to  a  chain  of 

three  vertices. 

Input  The  multilist  structure.  The  vertices  of  the 

polygon,  and  their  pointers. 

Output  The  updated  multilist  structure.  Updated 

degrees  and  cardinalities.  A  chain  in  stack 
CHAINC*).  Updated  value  of  M. 

2.5.19.  SUBROUTINE  REDUC4  (V  ) 

Description  This  routine  is  called  from  subroutine 

PY3467(V),  when  a  polygon  of  type  4  is  found. 
It  reduces  the  polygon  found  to  a  chain  of 

four  vertices. 


Input  The  multilist  structure.  The  vertices  of  the 

polygon,  and  their  pointers. 

Output  The  updated  multilist  structure.  Updated 

degrees  and  cardinalities.  A  chain  in  stack 
CHAIN(#).  Updated  value  of  M. 

3.5.20.  SUBROUTINE  REDUC6(V) 

Description  This  routine  is  called  from  subroutine 

PY3467CV),  when  a  polygon  of  type  6  is  found. 
It  reduces  the  polygon  found  to  a  chain  of 
four  vertices. 

Input  The  multilist  structure.  The  vertices  of  the 

polygon,  and  their  pointers. 

Output  The  updated  multilist  structure.  Updated 

degrees  and  cardinalities.  A  chain  in  stack 
CHAIN(*).  Updated  value  of  M. 

3.5.21.  SUBROUTINE  REDUC7(V) 

Description  This  routine  is  called  from  subroutine 

PY3467(V),  when  a  polygon  of  type  7  is  found. 
It  reduces  the  polygon  found  to  a  chain  of 
four  vertices. 

Input  The  multilist  structure.  The  vertices  of  the 

polygon,  and  their  pointers. 

Output  The  updated  multilist  structure.  Updated 

degrees  and  cardinalities.  A  chain  in  stack 
CHAIN(*).  Updated  value  of  M. 

3.5.22.  SUBROUTINE  P0LY2P(V) 

Description  This  routine  checks  if  vertex  V,  and  the  pair 

of  adjacent  vertices  specified  in  subroutine 
FINDNK(V)  form  a  polygon  of  type  2,  when  V  is 
not  in  set  K.  If  a  polygon  of  type  2  is 
found,  this  routine  calls  routine  REDC2P(V), 
which  performs  the  polygon-to-chain  reduction. 

Input  The  multilist  structure.  Vertex  V  and  a  pair 

of  vertices  adjacent  to  V,  along  with  their 
corresponding  degrees.  Pointers  to  the  pair 
of  adjacent  vertices,  in  V’s  list. 

Output  Indicator  FOUND  is  activated  if  a  polygon  of 

type  2  is  found. 


2.5.23.  SUBROUTINE  REDC2P  (V  ) 

Description  This  routine  is  called  from  subroutine 
P0LY2P(V),  when  a  polygon  of  type  2  is  found 
emanating  from  vertex  V,  when  V  is  not  in  set 
K.  It  reduces  the  polygon  found  to  a  chain  of 
three  vertices. 

Input  The  multilist  structure.  The  vertices  of  the 

polygon,  and  their  pointers. 

Output  The  updated  multilist  structure.  Updated 

degrees  and  cardinalities.  A  chain  in  stack 
CHAIN(*).  Updated  value  of  M. 

1*5.24.  SUBROUTINE  P0LY15(V) 

Description  This  routine  checks  if  vertex  V,  and  the  pair 

of  adjacent  vertices  specified  in  subroutine 
FINDNK(V)  form  a  polygon  of  type  1  or  5.  If  a 
polygon  of  type  1  or  5  is  found,  this  routine 
calls  routine  REDUCI(V)  or  REDUC5(V),  respec¬ 
tively,  which  performs  the  polygon-to-chain 
reduction . 

Input  The  multilist  structure.  Vertex  V  and  a  pair 

of  vertices  adjacent  to  V.  Pointers  to  the 
pair  of  adjacent  vertices,  in  V's  list. 

Output  Indicator  FOUND  is  activated  if  a  polygon  of 

type  1  or  5  is  found.  All  vertices  in  the 
polygon  found,  together  with  pointers  indicat¬ 
ing  their  addresses,  are  passed  to  the 
polygon-to-chain  reduction  routines. 

1*5.25.  SUBROUTINE  REDUC1 (V) 

Description  This  routine  is  called  from  subroutine 

P0LY15(V),  when  a  polygon  of  type  1  is  found. 
It  reduces  the  polygon  found  to  a  chain  of 
three  vertices. 

Input  The  multilist  structure.  The  vertices  of  the 

polygon,  and  their  pointers. 

Output  The  updated  multilist  structure.  Updated 

degrees  and  cardinalities.  A  chain  in  stack 
CHAINS).  Updated  value  of  M. 

1*5.26.  SUBROUTINE  R£DUC5(V) 

Description  This  routine  is  called  from  subroutine 

P0LY15(V),  when  a  polygon  of  type  5  is  found. 
It  reduces  the  polygon  found  to  a  chain  of  two 
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or  four  vertices,  depending  on  the  cardinality 
of  set  K. 

Input  The  multilist  structure.  The  vertices  of  the 

polygon,  and  their  pointers. 

Output  The  updated  multilist  structure.  Updated 

degrees  and  cardinalities.  A  chain  in  stack 
CHAIN ( * ) .  Updated  value  of  M. 

3.5.27.  SUBROUTINE  CHNRED 

Description  This  routine  performs  a  complete  reduction  on 
the  chain  resulting  from  a  polygon-to-chain 
reduction,  when  this  reduction  has  not  led  to 
a  cycle.  Vertices  V  and  W  are  the  chains  end¬ 
points.  The  routine  analyzes  three  cases:  a) 
deg(V)  >  2  and  deg(W)  s  2;  b)  deg(V)  =  2  and 

deg(W)  >  2;  and  c)  deg(V)  >  2  and  deg(W)  >  2. 

When  needed,  series  and  degree  2  reductions 
are  performed  and  the  T  list  is  updated. 

Input  The  multilist  structure.  The  chain  in  th 

CHAIN ( * )  array.  The  degrees  of  the  vertices. 
The  ONLIST ( * )  array.  The  T  list. 

Output  The  updated  multilist  structure.  The  updated 

T  list. 

1*5.28.  SUBROUTINE  OUTGRF 

Description  This  routine  optionally  prints  out  the  network 
when  it  has  not  been  totally  reduced  to  a 
two-vertices-one-edge  graph.  A  data  file  is 
always  written  containing  M  and  the  reduced 
network. 

Input  The  multilist  structure.  The  current  value  of 

M. 

Output  A  report  is  generated,  if  so  desired,  contain¬ 

ing  the  edges  of  the  reduced  network,  along 
with  each  edge  reliability.  In  this  report, 
the  current  value  of  M  is  also  listed.  A  data 
file  with  M  and  the  reduced  network. 

^.5.29.  SUBROUTINE  HEAD1 

Description  This  subroutine  prints  the  page  header  of  the 
reduced  network  output. 


Input 


Page  number,  date,  year,  and  hour. 


Output 


Page  header  with  incremented  page  number,  and 
date,  year,  and  hour  of  execution. 


3.5.30.  SUBROUTINE  OUTREL 

Description  This  routine  prints  out  the  network  when  it 
has  been  totally  reduced  to  a  two-ver tices- 
one-edge  graph. 


Input 


Output 


The  multilist  structure.  The  current  value  of 
M. 

A  report  is  generated,  containing  the 
network's  reliability. 


3.5.31_.  SUBROUTINE  HEAD2 

Description  This  subroutine  prints  the  page  header  of  the 
intermediate  and  final  results  output. 


Input 

Output 


Page  number,  date,  year,  and  hour. 

Page  header  with  incremented  page  number,  and 
date,  year,  and  hour  of  execution. 


.6.  Deb 


Feature  for  Program  Maintenance 


A  desirable  feature  in  any  code,  is  the  possibility  of  hav¬ 
ing  intermediate  printouts  of  important  variables  and 
arrays,  so  that  maintenance  can  be  simplified.  In  this 
code,  most  SUBROUTINES  have  a  LOGICAL  variable  named  TEST 
which,  when  "turned  on",  i.e.  when  made  equal  to  .TRUE., 
will  activate  such  a  debugging  feature. 
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4 .  User  * s  Guide 


In  this  section  of  the  report,  a  user's  manual  is  presented. 
First,  we  present  a  guide  for  using  PolyChain  in  an  applica¬ 
tion.  The  input  file  and  outputs  are  then  described. 
Finally,  two  test  problems  are  run  on  the  code  to  illustrate 
outputs . 

4_._1_.  Executing  PolyChain 

Here,  we  comment  on  how  the  program  is  executed.  The  algo¬ 
rithm  requires  that  the  network  be  nonsepar able .  An  algo¬ 
rithm  that  tests  for  this  condition  using  Depth  First  Search 
is  found  in  C A 7 4  3 .  To  run  the  code,  first  the  dimension 
parameters  and  the  COMMON  blocks  must  be  adjusted.  Then,  an 
input  data  file  must  be  prepared.  Both  of  these  topics  are 
presented  below. 

_4.J_.J_.  Dimension  Parameters 

The  first  step  in  running  PolyChain,  is  adjusting  the  dimen¬ 
sion  parameters  SUBROUTINE  INILST. 

These  parameters  are:  MAXVRT  -  maximum  number  of  vertices 

MAXEDG  -  maximum  number  of  edges 

If  the  dimension  of  the  problem  being  run  is  smaller  than 
the  dimensions  already  specified  in  this  subroutine,  no 
adjustment  is  needed.  But  if  the  network's  dimensions 
exceed  what  has  been  specified,  the  values  of  MAXEDG  and 
MAXVRT  must  be  altered.  Note  that  parameter  MAXCHN  need  not 
be  modified. 

After  adjusting  the  dimension  parameters,  all  COMMON  blocks 
containing  arrays  must  be  changed  accordingly.  These  COMMON 
blocks  are  adjusted  with  the  dimensions  specified  below. 


C0MM0N/BLK01 /  DEG (MAXVRT  ) 

COMMON/BLK02/  PTRADJ (MAXVRT ) , ADJVRT (2«MAXEDG ) , AVSADJ 

COMMON /BLK21 /  LNKDWN(2*MAXEDG),LNKUP(2*MAXEDG),LNKEDG(2*MAXEDG) 

COMMON /BLK03/  EDGPRB ( MAXEDG ), EDGNUM ( MAXEDG  ) 

COMMON /BLK04/  TLIST (MAXVRT ), CHAIN (MAXCHN ) ,TTOP,TOP, 

ONLIST (MAXVRT  ) 


4_.  J.  2  .  Input  Files 

Inputing  data  is  very  simple  in  PolyChain.  All  input  is 
format  free,  i.e.,  data  is  not  restricted  to  specific 
columns  of  the  input  line.  No  flag  is  needed  to  indicate 
end-of-f ile . 
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X 


The  first  line  of  the  input  file  contains  the  system  output 
options.  There  are  two  entries  on  this  line:  ECHOIN  and 
ECHOUT.  Set: 


ECHOIN 


{1  -  if  a  report  of  the  input  network  is  desired 
0  -  otherwise 


tv 


ECHOUT 


r  1  -  if  a  report  of  the  reduced  series-parallel 
complex  network  is  desired 

w 0  -  otherwise 


Setting  ECHOUT  to  "0"  will  not  hinder  PolyChain  from  writing 
a  data  file  with  M  and  the  reduced  network,  when  the  input 
network  is  series-parallel  complex. 


The  edges  are  entered  one  by  one  on  each  line.  To  specify 
an  edge,  enter  both  vertices  of  the  edge  followed  by  the 
edge's  reliability.  If  a  vertex  is  distinguished,  that  is, 
belongs  to  set  K,  it  should  be  preceded  by  a  minus  sign. 
Below,  we  illustrate  an  input  file,  for  the  network  i^  fig¬ 
ure  III. 


.v 

.v 
■  v 


Figure  III  -  Network  for  Example  Input  File 


k. 


For  this  network,  the  input  file  is  given  below.  Here 
wish  no  input  network  report,  but  require  a  report  of 
reduced  series-parallel  complex  network. 


0  1 

-1  2.8 

-1  3.3 

2  3  .2 


ft 


we 

the 
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2  -4  .7 

3  -4  .9 


4_.1_*1*  Program  Outputs 

In  this  section  PolyChain  is  tested  with  two  different  test 
problems  illustrating  the  program's  outputs. 

il*  K  Problem  One 

Consider  a  modification  of  the  ARPA  Computer  Network  in  fig¬ 
ure  IV. 


Figure  IV  -  Modified  ARPA  Computer  Network 


Consider  edge  reliabilities  for  this  network  varying  between 
0.1  and  0.9.  Both  output  options  are  set  to  "1".  The  input 
file  for  this  network  is  given  below. 

1  1 

-12.1 
-1  3  .2 
2  3  .3 

2  4  .4 

2  6  .5 

3  5  .6 

4  5  .7 

5  8.8 

5  20  .9 

6  11  .1 
6  20  .2 
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6  7  .3 

7  10  .4 

8  9  .5 

10  15  .7 

11  12  .8 

1213  .9 

1314  .1 

13  -21  .2 

14  15  .3 

9  16  .4 

16  17  .5 

1718  .6 

1819  .7 

19  -21  .8 


The  modified  ARPA  network  is  series-parallel  reducible. 
PolyChain  was  run  with  the  above  input  data.  For  this  type 
of  network  the  code  generates  a  two  part  report.  The  first 
part  describes  the  input  network,  edge  by  edge,  with  edge 
numbers,  vertex  numbers  and  types  (K=K-vertex,  nK=non  K- 
vertex),  and  edge  reliability.  This  section  also  summarizes 
the  input  network  data  and  core  usage.  The  second  part  of 
the  report  indicates  that  the  network  is  series-parallel 
reducible,  contains  the  K-terminal  network  reliability,  sum¬ 
marizes  the  reductions  performed  by  the  algorithm  and  indi¬ 
cates  the  CPU  time  (excluding  I/O).  The  three  page  report 
generated  by  PolyChain  for  the  above  file  follows. 
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Input  Network 
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— 

- - 
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9 
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5 

1 

1 
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20 

1 

1 

nK 
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nK 
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nK 
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Page  1 


Tue  Jul  26  1983 
12:59:52 


Reliability 


— — —————————  , 

. 10000000d+00  1 

I 


. 20000000d+00 


.30000000d+00 


. 40000000d+0  0 


.50000000d+00 
. 60000000d+00  | 
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Page  2 


Date  :  Tue  Jul  26  1983 
Time  :  12:59:52 

Input  Network 


Edge  ii  Vertex  i  Type  ii  Vertex  I  Type  i!  Reliability 

- ++ - + - +  + - + - +  + - 

- +  + - + - ++ - + - ♦+ - 

22  SI  16  nK  5  17  i  nK  i  |  . 50000000d+00 

—  - — ++- - - — + - ++ - ♦ - - 

23  ! !  17  I  nK  ii  18  i  nK  i  i  .  60000000d  +  00 

- ++ - — - + - ++ - + - ++ - 

24  | !  18  i  nK  ii  19  i  nK  i i . 70000000d+00 

- ++• - + - ++ - + - ++ - 

25  i  1  19  i  nK  ii  21  i  K  i  i . 80000000d+00 


Summary  of  Input  Network  Data 

Number  of  Vertices .  21 

Number  of  Edges .  25 

Number  of  K-Vertices .  2 

Network  Density .  0.119 

Summary  of  Core  Usage 

Variable  Current  Usage  % 

Name  Value 


MAXEDG 

MAXVRT 


5000 

2000 


25 

21 


0.5 

1.0 


Network  Ser ies-Parallel  Reducible 


Network  Reliability .  0. 87028600d-02 


Reductions  Performed 


Series .  19 

Degree  2 .  0 

Type  1 .  3 

Type  2 .  0 

Type  3 .  0 

Type  4 .  0 

Type  5 .  0 

Type  6 .  0 

Type  7 .  0 

Type  8 .  1 


Solution  Time  s  0.02  Secs. 


4.1_.3.2.  Problem  Two 


19  -21  .8 


For  series-parallel  complex  networks,  PolyChain  also  gen¬ 
erates  a  two  part  report  together  with  an  output  file.  The 
first  section  of  the  report  is  identical  to  the  series- 
parallel  reducible  case,  that  is,  contains  the  input  net¬ 
work.  The  second  part  indicates  that  the  network  is 
series-parallel  complex,  contains  the  reduced  network,  the 
updated  value  of  M,  and  summarizes  the  reductions  performed, 
including  percentage  reductions  in  the  network's  dimensions. 
CPU  time  (excluding  I/O)  is  also  indicated.  The  four  page 
report  generated  by  the  code  for  the  ARPA  computer  network 
is  given  next. 
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Summary  of  Input  Network  Data 


Number  of  Vertices . 21 

Number  of  Edges .  26 

Number  of  K-Vertices .  2 

Network  Density . . .  0.124 


Summary  of  Core  Usage 


Variable  Current  Usage  % 

Name  Value 


HAXEDG 

MAXVRT 


5000 

2000 


26 

21 


0.5 

1.0 
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Date  :  Tue  Jul  26  1983 
Time  :  12:58:26 


Updated  value  of  M  =  0 . 21 67 9720d+00 


Reductions  Performed 


Series .  15 

Degree  2 .  0 

Type  1 .  3 

Type  2 .  0 

Type  3 .  0 

Type  4 .  0 

Type  5 .  0 

Type  6 .  0 

Type  7 .  0 

Type  8 .  0 


Original  Reduced  %  Reduction 

Network  Network 


Edges .  26  8  69.2 

Vertices .  21  6  71.4 

K-Vertices .  2  2  o* 


Solution  Time 


0.02  Secs 
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For  the  series-parallel  complex  case,  PolyChain  also  writes 
a  file  with  the  updated  value  of  M  and  the  reduced  network. 
We  next  present  the  file,  named  "polychain  .out" ,  generated 
in  this  example.  This  file  can  be  used  as  input  to  a  fac¬ 
toring  algorithm  program. 


0.2167  9720d+00 


-1 

9 

0.31 42281 2d+00 

-1 

6 

0.5381 7522d+00 

6 

13 

0. 18000000d-01 

6 

10 

0 . 40000000d-0 1 

9 

-21 

0. 1 1760000d+00 

9 

10 

0.600000006+00 

10 

13 

0. 64000000d-01 

13 

-21 

0. 30000000d+00 

The  reduced  network  obtained  by  PolyChain  for  the  ARPA  com¬ 
puter  network  is  given  in  figure  VI. 


Figure  VI 


Reduced  Network 
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5.  Large-Scale  Networks 


Next,  we  test  PolyChain  with  a  number  of  large-scale  net¬ 
works.  These  networks  help  us  understand  PolyChain's 
behavior  when  treating  large  networks  and  enable  the  estima¬ 
tion  of  CPU  times  for  even  larger  problems.  A  random  net¬ 
work  generator  was  used  to  provide  the  series-parallel  com¬ 
plex  networks,  while  the  series-parallel  reducible  networks 
were  built  aggregating  a  number  of  150  vertex,  208  edge  net¬ 
works.  Problems  were  run  on  the  VAX  11/750  of  the 
Etcheverry  Hall  VAX/UNIX  CAE  Laboratory,  at  Berkeley.  The 
code  was  compiled  on  the  UNIX  f77  compiler  using  the  -0 
optimizing  option.  The  code  was  run  under  the  UNIX  operat¬ 
ing  system.  CPU  times  were  measured  through  the  "dtime" 
system  routine.  Table  I  contains  a  summary  of  the  networks 
tested . 
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K-Ver tices 

Density 
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63 
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No 
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Deg2 
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16 

49 
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1  0 

m 

2 

35 
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32 

99 

12 

m 

20 

m 

■ 

71 
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48 
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18 

m 

30 

■ 

6 

107 
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96 
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36 

_ 

m 

60 

n 

12 

215 

98  18  1  1 


216  50 


330  36 


583  101 


Mean 

CPU 

Time 


0, 

.27s 

0, 

,57s 

1  , 

,09s 

1  , 

.  60s 

3. 

,23s 

0, 

.30s 

0. 

.  60s 

1  , 

.14s 

1  , 

,  66s 

3. 

.58s 

Table  II  -  Test  Results 


Each  test  problem  was  run  10  times.  The  problems,  however, 
do  not  cover  a  very  wide  spectrum  of  networks  and  therefore 
some  bias  may  exist  in  our  conclusions. 

The  tests  agree  with  the  classification  of  the  algorithm  as 
being  O(iEi).  The  network  size  (!E!)  processed  per  CPU 
second  varied  from  733.33  to  784.28  edges/sec  for  the 
series-parallel  reducible  networks  tested  and  from  694.08  to 
766.26  edges/sec  on  the  series-parallel  complex  networks. 
The  mean  network  size  <!Ei)  processed  per  CPU  second  was, 
respectively,  775.37  and  725.17  edges/sec  for  the  reducible 
and  complex  networks. 

These  first  results  show  that  PolyChain  is  feasible  for  com¬ 
puting  K-terminal  reliability  of  large-scale  series-parallel 
reducible  networks.  With  1  min.  CPU,  one  should  be  able  to 
compute  the  K-terminal  reliability  of  a  series-parallel 
reducible  network  upwards  of  35,000  edges. 


6.  Conclusions  and  Recommendations 


6.1.  Conclusions 


This  report  discussed  the  design  and  implementation  of 
PolyChain,  a  portable  FORTRAN  code  for  evaluating  undirected 
network  reliability  via  polygon-to-chain  reductions,  [S82], 
[W82].  The  code’s  implementation  facilitates  further  exten¬ 
sions  and  enhancements.  It  uses  a  multilist  data  structure 
representation  of  the  network,  which  enables  good  core 
management  and  efficient  network  manipulation.  Input  is 
simple  and  format  free.  Input  data  is  tested  for  con¬ 
sistency  and  execution  is  terminated,  when  inconsistent  data 
is  detected.  Outputs  are  detailed  and  a  data  file,  contain¬ 
ing  the  reduced  network,  is  written,  when  the  network  is 
series-parallel  complex. 

The  code  was  tested  on  large-scale  networks  and  the  results 
are  encouraging. 

A  copy  of  the  program’s  source  code  is  available  through  the 
author . 


6.2.  AMCOModations 


Further  testing  is  still  needed  to  ensure  the  code’s 
correctness.  A  random  network  generator  is  also  required, 
so  that  large-scale  networks  can  be  further  tested  on 
PolyChain.  This  generator  should  have  the  capability  of 
generating  both  series-parallel  reducible  and  complex  net¬ 
works  . 

The  algorithm  has  a  constraint  on  the  topology  of  the  input 
network.  The  network  must  be  nonseparable .  In  the  present 
version,  PolyChain  does  not  test  for  this  requirement. 
Since  these  networks  may  occur  in  practice,  a  routine  to 
check  for  this  condition  is  desirable.  Aho ,  Hopcroft,  and 
Ullman,  [A7M],  present  a  Depth  First  Search  based  algorithm 
for  determining  all  nonseparable  components  of  a  network. 

To  insure  the  evaluation  of  the  K-terminal  network  reliabil¬ 
ity  for  any  network,  this  code  should  encorporate  a  factor¬ 
ing  algorithm.  In  this  manner,  the  program  could  oscillate 
between  both  algorithms  until  a  final  result  is  attained. 
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